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The dominant theoretical uncertainties in both, the 
anomalous magnetic moment of the muon and the value of 
the electromagnetic coupling at the Z scale, Mz, arise from 
their hadronic contributions. Since these will ultimately dom- 
inate the experimental errors, we study the correlation be- 
tween them, as well as with other fundamental parameters. 
To this end we present analytical formulas for the QCD contri- 
bution from higher energies and from heavy quarks. Including 
these correlations affects the Higgs boson mass extracted from 
precision data. 

PACS numbers: 13.40. Em, 14.60.Ef, 12.20.Ds, 12.38.Bx. 



The magnetic moment of the electron, g e , provides the 
best determination of the fine structure constant, but 
is currently not measured precise enough to give a sen- 
sitive probe of electroweak physics which is suppressed 
by a factor {a/ir)m 2 e /M^ ~ 0(l(r 13 ). Meaningful (but 
still weak) bounds on new physics contributions can be 
set only if they enter without loop suppression. On the 
other hand, electroweak contributions to the anomalous 
magnetic moment of the muon Q], = (py, — 2)/2, are 



and where 



enhanced by a factor m^/m^ ~ 4 x 10 4 , which renders 
them sizable enough to be detectable at the ongoing E821 
experiment at the AGS at BNL. (a T has not yet been 
observed experimentally). E821 already reduced the ex- 
perimental error to ±1.6 x 10~ 9 ||. The anticipated final 
error of about ±0.4 x 10~ 9 will mean a factor of 20 im- 
provement relative to previous results j3|. a M therefore 
provides a good laboratory to test the Standard Model 
(SM) and probe theories beyond it 0. For example, 
scenarios of low-energy supersymmetry with large tan j3 
and moderately light superparticle masses can give large 
contributions to a M ||. 

Unfortunately, the interpretation of is compromised 
by a large theoretical uncertainty introduced by hadronic 
effects. Two and three loop vacuum diagrams containing 
light quark loops cannot be calculated reliably in per- 
turbative QCD (PQCD). Instead they are obtained by 
computing dispersion integrals over measured (at low en- 
ergies) and theoretical (at higher energies) hadronic cross 
sections. At two loop M, 



a M (had; 2-loop) = 
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where R(s) is the cross section of e + e — ► hadrons, nor- 
malized to the tree level cross section of e + e~ — > 
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The uncertainty introduced by this procedure is signifi- 
cantly larger than the anticipated experimental one. For 
example, Ref. quotes an error of ±0.67 x 10~ 9 , while 
other evaluations [p|-|l(| give larger uncertainties. An 
analogous uncertainty occurs in the QED coupling con- 
stant, a(n), (MS quantities will be marked by a caret), 
preventing its precise theoretical computation from the 
fine structure constant, a, for // > 2771^0. Knowledge of 
a(Mz) is indispensable for the extraction of the Higgs 
boson mass, Mh, from the mass of the W boson, Mw, 
and the weak mixing angle, Syy- Again one must rely on 
a dispersion integral to be taken over R(s), but using a 
different kernel function. As a result, these hadronic un- 
certainties are strongly correlated with each other, and 
also with other fundamental SM parameters, such as the 
strong coupling constant and the heavy quark masses. 

To address these correlations we obtain compact ana- 
lytical expressions for a M (had) wherever possible, i.e. for 
the charm and bottom quarks, as well as for the high en- 
ergy integration region (above 1.8 GeV) for light quarks 
where PQCD appears to be applicable. A similar idea has 
been pursued in Ref. [jllj, for the calculation of a(Mz). 

To lowest order in QED, the anomalous magnetic mo- 
ments are lepton universal and mass independent |12| , 
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To two-loop order one has, 

2-loop = nl i , n e , n r , 



^uds 



^had 
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where the lepton contributions |13| amount to a^J 



butions. To proceed, we used the expansion, 



= 4132.18 x 10 . a u and a c u are heavy quark contri- 
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K(s) = 1 + fci(Vi)-^ + fc(>/S)-£ + 0(-£) 
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In Eq. (0) the integration is over the imaginary part of 
the photon polarization function which is related to R 
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by R(s) — 127rImII((7 2 + ie) and is taken along the real 
axis from the quarkonium threshold to infinity. Ana- 
lytic continuation allows us to integrate instead over the 
full function, n(s), along a circle of radius s — rh 2 (rh 2 ) 
counterclockwise around the origin. This avoids a com- 
plicated integration over resonances, i.e. a region where 
perturbative QCD is not applicable, at the expense of 
introducing a new uncertainty from the imperfect knowl- 
edge of the MS quark mass (see the discussion in the next 
paragraph). With II(s) known up to 0(a 2 s ) |l4| ], we find 
for a quark of charge Q q , 
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The remaining terms in Eq. (Q) are due to u, d, and s 
quark effects, which we separated into the contributions 
from V« > Mo = 1-8 GeV {af s ) and V« < Mo {a^)- 
a l ^ ds can be written as an expansion in to 2 /^ 2 ,. For the 
leading contribution we find, 
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where d 2 = 299/24 - 9C(3) and d 3 = 58057/288 - 
779C(3)/4 + 75C(5)/2 are coefficients of the Adler D- 
function, and 
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which we compute with a 4-loop renormalization group 
improvement. For a representative value a s (^o)/ 71 ' = 0.1, 
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where a s — a s (m q )/ir, ni = 3 for charm, and ni = 4 for 
bottom. Higher orders in to 2 /m 2 can safely be neglected. 
Terms of 0(a 2 TO^) can also be dropped unless they are 
logarithmically enhanced. Non-perturbative effects in 
the operator product expansion were also computed and 
found to be negligible. Using (here and in the following) 
a s (M z ) = 0.120, m c = 1.31 GeV, and rh b = 4.24 GeV, 
we find al = 1.39 x 10" 9 and al = 0.03 x 10~ 9 . 

In contrast to the numerical integration over reso- 
nances, Eq. (^) is a simple and transparent representation 
of the heavy quark contribution. More importantly, the 
uncertainty implied by Eq. (Q) is smaller. To see this, 
notice first the excellent behavior of its a s expansion im- 
plying a small truncation error which we evaluated to be 
±0.05 x 10~ 9 . The dominant uncertainty is induced via 
the MS charm mass, but it is necessarily smaller than 
in the conventional approach regardless of the available 
data in the resonance region: if one compares the two 
treatments — integration over the real axis heavily re- 
lying on experimental results and integration over the 
circle contour relying only on PQCD — one can deter- 
mine the heavy quark masses by demanding consistency. 
This effectively amounts to deriving a specific QCD sum 
rule. The point we are making here, is that this is only 
one of a large number of possible sum rules, but not the 
one which uses the available information most efficiently. 
The charm mass extracted from the most efficient sum 
rule fl5| (and other precise quark mass determinations) 
can then be used for Eq. ([7). An unweighted average 
of various determinations yields |ll[] Am c — ±0.07 GeV 
and Aaf. = ±0.16 x 10~ 9 . In the near future these un- 
certainties are likely to reduce to about ±0.04 GeV and 
±0.09xl0 -9 , respectively. Similarly, Am b = ±0.11 GeV, 
but the induced error in a M is negligible. 



Bi = 7.069 x 10" 
B-i = 2.562 x 10" 



B 2 = 4.514 x 10" d , 

B 4 = 1.243 x 10~ 5 . (10) 



Notice, that B4 is small enough, that even with the fourth 
order coefficient, d±, unknown this treatment keeps the 
truncation error at a negligible level. (For an estimate 
of the uncertainty of di, see Ref. Denoting k[ = 

ki(no) — § and k' 2 = /^(/io) — 6, we find for the subleading 
contributions, 
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where in 0(a s ) we kept the small s quark mass effect 
(~ 3 x 10~ 12 ). G(x) arises from virtual charm quark 
effects inside a light quark loop (double bubble diagram) . 
Even below threshold it can be well approximated as an 
expansion in x despite the fact that m 2 < /1q, 
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G(x) also applies to b quarks, but this contribution can 



be safely neglected. We find a) 



uds 



4.38 x 10" y and 



the leading order resummation in Eq. (g) renders the 
truncation error negligible. 

Leading order non-perturbative contributions due to 
gluon and light quark condensates are given by, 
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where the condensates are of order Aq CD and — m|-/^, 
respectively. These terms are suppressed by two pow- 
ers of a s and they change et M by less than 10 -12 , an 
effect completely negligible. Effects from up and down 
quark condensates are suppressed by a further factor of 
m^/m 2 K i and quartic mass terms are tiny, as well. Un- 
certainties from non-perturbative effects not accounted 
for by the operator product expansion, which are due to 
the transition from the data region to PQCD at no, have 
been estimated in Ref. which quotes ±0.024 x 10~ 9 . 
We take the low energy contribution from Ref. M , 
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(63.43 ±0.60) x 10" 
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which includes a QCD sum rule improvement, PQCD 
down to a relatively low = 1-8 GeV, as well as ad- 
ditional information from r-decays. The quoted error 
is not uncontroversial and needs to be confirmed. 

Note, however, that inclusion of the r-decay data de- 
creases the difference between the SM prediction and the 
current experimental result Q . Our reason to use Ref. 0] 
is that it quantifies the correlation (69%) with the cor- 
responding result on Aahad = (56.53 ± 0.83) x 10~ 4 . 
We also account for the almost perfect anti-correlation 
with higher order hadronic uncertainties, as will be dis- 
cussed below. The total two-loop quark contribution is 



a b „ + a?, + a uds 



> ad = 69.23 x 10~ 9 . 
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This completes the discussion of the terms appearing 
in Eq. (|4|). The leading contribution at 0(a 3 ) is from the 
light-by-light diagram containing an electron loop JlTj , 
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where we displayed only the leading term in explicit 
form, but the terms suppressed by one or two powers 
of rrie/m^ are significant, as well. The contribution due 



to a' 



;lbl 



20.9479 ot f 7r = 262.54 x 10" 9 is almost 
four times larger than the entire hadronic contribution 
at 0(a 2 ). The corresponding contribution involving a 
muon loop is a£ ;lbl = 0.37f0 a 3 /n 3 = 4.65 x 10" 9 , and 
the r contributes [fL7| , 
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An evaluation of the hadronic light-by-light contribution 
yields p|, 



z had;ibi = (_o.792 ± 0.154) x 10" 



(16) 



which is consistent with the finding in Ref. |l9fl . No- 
tice, hadrons and leptons contribute with opposite signs. 
Since these results are based on model calculations, an 
independent confirmation would be desirable. 

The purely leptonic 3-loop vacuum polarization con- 
tribution is again dominated by electron loops and given 



by a^ v P Ql = 2.7294 a 3 /ir 3 = 34.21 x 10~ 9 . On the 
other hand, the contribution from two hadronic loops is 
suppressed by a factor m 4 t /16m 4 . ± and therefore small, 
a had;v P oi = (0.027 ± 0.001) x 10~ 9 §. Similar to our 
strategy at two loops, we separated the mixed leptonic- 
hadronic contribution into heavy quarks, light quarks 
(V* > Mo); an d light hadrons (y/s < /i ). The r-hadronic 
contribution is of the order of 10~ 12 and negligible. As 
for the other leptons, we use the kernel functions from 
Ref. pel and obtain for charm and bottom quarks, 
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which amounts to —0.05 x 10 9 and —0.002 x 10 9 , re- 
spectively. For the light quarks we obtain, 

a 3 1 + a s 



£—uds;vpol 





9 



/ 20359 103 



4^ 



509 



23, mi 
3 ^5 



(4 \ 576 



^ + ^7 ln ^+ (18) 



24 



72 



19 m£ 
6 mt 



2-1 - ii 2 



■In' 



m 2 m 2 
2 In ^ In — % 

Mo m » 



/'o 

= -0.15xl0~ 9 . 



We found the vacuum polarization contribution aris- 
ing from an electron in one loop and light hadrons in 
the other one, a ^- had;vpo1 = 0.97 x 10~ 9 , by constructing 
a simplified function R(s) which reproduces the results 
in Ref. §. a e-had ; v P oi is a i most completely (> 99.9%) 

correlated with a£ ad in Eq. @. On the other hand, 

a M-had;vpol = _ 18 q x 1Q -9 jg yery strong l y (~ -97%) 

anti-correlated with a^ ad . The small uncorrelated error 

a* 

contributions are clearly negligible, and the correlated 
ones can be added (subtracted) linearly, slightly reduc- 
ing the error in Eq. ( [l4| ) to ±0.59 x 10~ 9 . Including 
the uncertainty from the light-by-light contribution in 
Eq. ( [llf ) we obtain ±0.61 x 10~ 9 as the total hadronic 
error excluding parametric uncertainties. 

Taking the other theoretical uncertainties mentioned 
earlier to be 100% correlated with the corresponding un- 
certainties in Aahad (which is a fit parameter) , we obtain 
a residual correlation of 55% or ±0.34 x 10 -9 (the uncor- 
related error is ±0.51 x 10~ 9 ). 

The four-loop contribution |§l), a^ loop = 126.04a 4 /7r 4 
= 3.67 x 10~ 9 , and the five-loop estimate ||, a^~ loop = 
930a 5 /7r 5 = 0.06 x 10~ 9 , are also included. 

One-loop electroweak corrections due to W and Z bo- 
son loops are given by p3[, 
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Two-loop corrections to Eq. ( [19|) are significant due to 
large logarithmic contributions |24|-Hffi|], 
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The fermionic two-loop result, including some additional 
logarithms and sub-leading contributions, has been ob- 
tained in Ref. j2fj for small and large values of Mh- Wc 
constructed an interpolation formula for other values of 
Mh which reads in units of a(m /J )G_Frn^/(3v / 27r 3 ), 
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where x = m 2 t /M H and w = e ~ 3 - 31849x . Eq. @ re- 
produces the exact result |^5| for Mh — rh t . The 
bosonic two-loop corrections have been obtained only 
for large values of Mh and as an expansion in s\j pfjj. 
We take the leading non-logarithmic contribution in the 
Mh — » oo limit (a_2 in Ref. p6[ ) as the uncertainty 
induced by sub-leading bosonic two- loop effects. The 
leading logarithms to three-loop order have also been 
computed |27|] and included in our analysis. We obtain 



,EW _ 



(1.52 ±0.03) x 10" 



Summation of the various contributions gives (g^ — 2 — 
-)/2 = 4506.28 x 10~ 9 . Including our evaluation of a M 
into a global analysis of electroweak data yields, 

-0„ -2--) = (4506.35 ± 0.37 ± 0.51) x 10~ 9 , (22) 

2 7T 

where the first error includes Aahad and all other para- 
metric uncertainties, such as the ±0.028 error in a s {no) 
obtained from current global fits. The extracted Higgs 
mass from this fit is Mh — 88^33 GeV. The current ex- 
perimental world average, (4510.55 ± 1.51) x 10~ 9 (which 
differs by 2.6cr from Eq. (|22|)) needs to be included into 
the analysis, as well. This is achieved by combining the 
experimental error with the second error in Eq. (E2J) in 
quadrature (for a total of 1.59x 10~ 9 ). The global fit then 
yields 4506.52 ±0.36 (parametric) with a pull of 2.5. The 
crucial observation here is that now M H = 83lg[ GeV, 
i.e. both, the central value and the uncertainty of Mh 
decrease. The increase in the precision in Mh gained by 
properly correlating it to Aahad (and assuming the SM) 
is almost identical to that provided by the Myy measure- 
ment at the Tevatron Run I. 

To summarize, we carefully examined the error corre- 
lations between and other quantities entering elec- 
troweak tests, especially Aa(Mz), a s , and the quark 
masses. We derived new analytical results for the 
hadronic contributions and showed that the proper treat- 
ment discussed in this article has a significant effects on 
the extraction of Mh, which could (depending on the fu- 
ture experimental findings) become even more dramatic. 
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